Model 2 analysis
Benjamin Mercier
Last updated on 16 November, 2022
Model 2
Description of the model
This second model is a slight variant of the first, in which we assign a specific \(\alpha\) to each predator. In doing so, we can isolate difference between predator, to check if some predator seems to be different from the whole. \[\begin{align}
F_{ij}^{real} &= \alpha_{j} * B_i * \frac{B_j}{M_j}
\end{align}\]
This model was fit with a hierarchy implemented on the alpha parameter. A global alpha was estimated, with 118 respective unique alphas for each predators.
Summary table
| mean | se_mean | sd | 2.5% | 25% | 50% | 75% | 97.5% | n_eff | Rhat | |
|---|---|---|---|---|---|---|---|---|---|---|
| a_pop | -10.218115 | 0.0028599 | 0.3273803 | -10.854945 | -10.443121 | -10.218501 | -9.998045 | -9.582319 | 13104.332 | 0.9998744 |
| a_sd | 3.468383 | 0.0021624 | 0.2303778 | 3.055874 | 3.306278 | 3.453411 | 3.614110 | 3.964302 | 11350.477 | 0.9997395 |
| sigma | 1.682513 | 0.0004158 | 0.0406853 | 1.606718 | 1.654718 | 1.680976 | 1.709375 | 1.764458 | 9572.696 | 0.9997438 |
Exploring parameters posterior distribution
Respective predator alphas
One-one plot of simulation against data
